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Abstract 

Allozymes and mtDNA sequences were used to assess genetic variability in six hunted and three non- 
hunted populations of red foxes Vulpes vulpes in Central Europe and the Mediterranean area. Allo- 
zyme variability was very low (P = 0.0- 4.4; He = 0.0-0.006) and significantly smaller in protected pop¬ 
ulations than in hunted ones (P = 8.9-15.6; He = 0.010-0.044). As a trend, this result was confirmed by 
the distribution of mtDNA variation. Differences of genetic variability were independent from sample 
size. Most likely, foxes evolved under predation by larger carnivores, which probably have primed the 
turnover of territorial reproducers. In absence of predation, hunting could superficially mimic such an 
effect and favour a greater genetic variability and, thus, a better viability in the red fox. This may sug¬ 
gest an explanation to the failure of intensive hunting campaigns as a measure to reduce strongly the 
population size of red foxes. 
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Introduction 

The social organisation of the red fox Vulpes vulpes (L., 1758) includes territorial individ¬ 
uals, who are resident and thus responsible for reproduction, as well as itinerant individ¬ 
uals, who are subordinate and often subadult, who may reproduce only when they replace 
a territorial fox (Macdonald 1980 a, 1987). In a fox population genetic variability is most 
likely dependent on effective population size (Ne), on frequency of replacement events in 
territories by itinerant individuals, and on the relatedness of the replacing individual to 
the group. Replacement rate should increase in populations with a high mortality rate in 
the adult age class, induced by predation or hunting. Several proximate evolutionary de¬ 
terminants may also have an effect on reproduction of red fox populations, thus presum¬ 
ably on their genetic structure. These include effects of mortality rate on the proportion 
of reproducing vixens and on litter size, as well as dispersal “sinks” (e. g. Harris and 
Smith 1987; Zimen 1984; Artois 1989, for a review). 

Allozymes and DNA sequences are useful markers to study population genetics in 
mammals. Statistical indices inferred from allele frequencies at protein loci and the distri¬ 
bution of different mitochondrial DNA haplotypes can be used to quantify the extent of 
genetic variability and differentiation of natural populations. Estimates of genetic diver- 
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sity can be correlated with several parameters (ecology, climate, vegetation, population 
density, life cycle and style), which may influence the level of variability (Nevo et al. 
1984), thus providing useful suggestions for management (Moritz 1994). The organisation 
of mating systems (Apollonio and Hartl 1993) and other factors, e. g., replacement rate 
of territorial individuals where they are the only reproducers, may also influence genetic 
variation. In this study we examine the relationships between hunting and protection 
from hunting with genetic variability in red fox populations. 


Material and methods 

Standard allozyme and DNA analyses have been carried out on samples of liver tissue removed from 
foxes freshly killed by hunters in the course of occasional sampling/targeted control operations (i. e. 
rabies monitoring and livestock protection) and regular hunting (Tab. 1) from selected sites in Central 
Europe and the Mediterranean range (Fig. 1). Sampling in protected areas was limited because of 
ethical and legal constraints on the collection of specimens, although all sampled populations of red 
foxes were abundant. All samples came from adult (> 1 year old) individuals to decrease the probabil¬ 
ity of close relatedness between them. To our knowledge, no bottleneck has occurred in these popula¬ 
tions in the last few centuries. Liver samples were preserved at -80 °C prior to analyses and used to 
perform an electrophoretic screening of 45 presumptive loci (Grillitsch et al. 1992; Hartl and Ho; 
GER 1986). A subset of the same samples was used to extract DNA. The 5'-end of the mitochondrial 
Cytochrome b gene (Cyt b: 375 bp) was PCR-amplified and sequenced according to the procedure 
outlined in Simon et al. (1991). A detailed analysis of several statistical parameters derived from allele 
frequencies and from DNA sequences was obtained with the computer programs BIOSYS-1 (Swof- 
FORD and Selander 1989), REAP (McElroy et al. 1992) and MEGA (Kumar et al. 1993) and was 
presented in Frati et al. (1998). 


SP2 ITl IT3 AUl BU2 BUI 



Fig. 1. Map of collecting sites for the nine fox populations: Spain (SPl, SP2), Italy (ITl, IT2, IT3), Aus¬ 
tria (AUl), Bulgaria (BUI, BU2), Israel (ISl). 
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The nucleotide sequences reported in this study have been deposited in the EMBL, GenBank and 
DDBJ Nucleotide Sequence Databases under accession numbers Z80957-Z80984 and Z80987-Z80997 
(cf. Frati et al. 1998). 


Results and discussion 

Several indices of genetic variability, inferred from both allozymes and mtDNA, are 
shown in table 1. Red foxes in non-hunted populations showed a significantly lower level 
of allozyme variability than those in hunted populations (Mann-Whitney U-test, two- 
tailed: U = 0.0; ni = 3; n 2 = 6; p < 0.05, in both the proportion of polymorphic loci and het¬ 
erozygosity; Tab. 1). Our results might also be a consequence of sample sizes, e. g., only a 
small number of family groups may have been sampled in the unhunted sites. In fact, 
there was a significant correlation between sample size and allozyme variability (r = 0.68, 
df = 7, p < 0.05). We reduced all sample sizes to three specimens, i. e., the smallest sample 
(ITS), randomly chosen within each population, to avoid partly such a correlation. In 
spite of this reduction (repeated three times, independently), hunted populations re¬ 
mained much more variable than non-hunted ones at the same level of significance and 
parameters as above (Tab. 1). 

A relationship between hunting and the increase of genetic variability appeared to be 
confirmed by the two indices inferred from Cyt b sequence data (Tab. 1), except for the 
ITS population, which showed a higher level of variability than that observed with allo¬ 
zymes. Sequence divergence between haplotypes ranged from 0 to 2.67% and it was max¬ 
imum among specimens of the AUl population. 

Our data suggest that genetic variability in unhunted fox populations is almost absent, 
independently from sample size, and it is significantly lower than that in hunted populations 
where heterozygosity and proportion of polymorphic loci are within the range of those 
found amongst other wild canids (Hamilton and Kennedy 1986; Kennedy et al. 1991; Lor- 
ENZiNi and Fico 1995). Differences in polymorphism are found not only at the population 
level, but also at the individual level (Tab. 1, Ho). This should rule out the possibility that a 
lower variability in non-hunted populations is due to sampling of related individuals. 

Variability of fox behaviour and ecology is influenced by habitat structure and rich¬ 
ness (e. g. Artois 1989; Lucherini and Lovari 1996; Macdonald 1980 a). This includes 
aspects of the reproductive system and of spatial organisation (territory size, dispersal 
patterns), which in turn may influence genetic structure. The most likely factor underlying 
an unusually low genetic variability in a population is the local scarcity of assortative mat¬ 
ing (Crow and Kimura 1970). Non-territorial foxes occur in non-hunted populations (e. g. 
Cavallini and Lovari 1994; Lovari et al. 1994; Lucherini et al. 1995), but apparently they 
fail to reproduce. Therefore, the nomadic fraction of a fox population may not contribute 
to reproduction, unless territories become vacant. This effect of protection is perplexing 
as natural, unhunted populations should be able to maintain enough genetic variability 
through ecological and behavioural mechanisms. 

Polygynous mating systems decrease genetic variability among ungulates, because only 
a few males contribute to reproduction (Apollonio and Hartl 1993). If foxes were poly¬ 
gynous, i. e. one male would reproduce with several females, this could favour the loss of 
genetic variability. Up to six vixens are known to share their territory with one dog fox, 
especially at the comparatively high density found in protected areas, with no predators, 
but normally only the alpha female has cubs (Macdonald 1979; Schantz 1981; Zabel 
and Taggart 1989). Therefore, little or no difference in allozyme variability can be ex¬ 
pected between polygynous and monogamous fox populations. 

There may be another explanation to the low level of genetic variability in protected 
fox populations. The general rule amongst canids is that “larger species dominate - or 


Table 1. Estimates of genetic variability in nine populations of red foxes. Abbreviations as in fig. 1. (Compiled from Frati et al. 1998). ss: sample size, P: Per¬ 
centage of polymorphic loci (99% criterion); He: Expected heterozygosity according to the Hardy-Weinberg equilibrium, Ho: Observed heterozygosity. 
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even kill - smaller ones” (Macdonald 1992: 90). This has been reported for, e. g., wolves 
and coyotes (Carbyn 1982; Fuller and Keith 1981), coyotes and red foxes (Dekker 
1983; Major and Sherburne 1987; Voigt and Earle 1983; but see also Gese et al. 1996), 
red and arctic foxes (Bailey 1992; Frafjord et al. 1989; Hersteinsson and Macdonald 
1992; Rudzinski et al. 1982), coyotes and kit foxes (Ralls and White 1995; White et al. 
1994). Most likely, if the red fox evolved under the pressure of a community of larger pre¬ 
dators such as other canids as well as felids, e. g., cougar and leopard (Cavallini 1996), 
this may have strongly influenced the development of its great biological flexibility and 
resilience (Nowak 1991). The only larger canid in southern Europe is the wolf, which is 
known to feed on foxes on occasion (Matteucci 1992; Meriggi et al. 1991). Presently, 
predation pressure by the wolf on foxes is likely to be small in most of Europe, as the 
wolf is an endangered species while the red fox is one of the most successful mammals in 
terms of numbers and range width. The Iberian lynx Lynx pardinus occurs in one (SPl) 
of the protected areas of our sample, but it is rare and unlikely to prey on red foxes be¬ 
cause of its small size (Rau et al. 1985). No wolves - or any other larger predator - have 
occurred in our other protected sampling areas for at least a few decades. The lack of na¬ 
tural predators may result in relatively stagnant fox populations, where the same indivi¬ 
dual could successfully defend a territory for several years. This may prevent a frequent 
turnover of reproducers over several years, eventually leading to a decrease of genetic 
variability. In polygynous groups of red foxes, vixens are closely related, often mothers 
(the alpha female) and daughters (Macdonald 1979; Schantz 1981). If the alpha female 
dies, one of the “helpers” will replace her as reproducer (Macdonald 1987). Thus, some 
inbreeding can be expected. The occurrence of protected areas, devoid of larger preda¬ 
tors, is a very recent event in fox evolution. No equally effective, alternative process to 
predation may have yet been developed in red foxes to favour genetic variability. 

Hunting may increase genetic variability in fox populations by partly mimicking the 
effects of natural predation with which foxes have probably evolved. A high genetic 
variability in a population is believed to enhance its viability (Allendorf and Leary 
1986; Lacy 1997, for a review). If our hypothesis is correct, this may suggest one impor¬ 
tant reason why persecution campaigns or intensive hunting have never succeeded to era¬ 
dicate foxes from large areas, e. g., for rabies control in Europe (Macdonald 1980 b, for a 
review). In this connection, the ecological adaptability of a generalist species, such as the 
red fox, should also not be underestimated. 

On the other hand, hunting alone is unlikely to be a valid replacement for natural pre¬ 
dation. Long-lasting hunting practices can alter the genetic structure of mammalian popu¬ 
lations, often acting in an anti-Darwinian sense (Hartl et al. 1991; Hartl et al. 1995; Ry- 
MAN et al. 1981; Scribner et al. 1985), thus being poor substitutes of natural predation. 
We wish to stress that our results need confirmation from larger samples and from obser¬ 
vational data of territory turnover in areas where a community of predators of foxes does 
and does not occur. If hunting is in principle an acceptable activity, there is no reason not 
to admit it for red foxes. Conversed, our data must not be used to justify intensive hunt¬ 
ing or persecution campaigns, which are quite a different, controversial issue on biological 
and ethical grounds, beyond the scope of this paper. 
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Zusammenfassung 

Trdgt die jagdliche Schonung zur genetischen Einformigkeit beim Rotfuchs bei? 

Mittels Allozymelektrophorese und mtDNA-Sequenzierungen wurde die genetische Variabilitat von 
sechs bejagten und drei nicht bejagten Populationen des Rotfuchses Vulpes vulpes in Mitteleuropa 
und in mediterranen Gebieten untersucht. Die genetische Variabilitat war deutlich geringer in den 
nicht bejagten (P = 0,0-4,4; He = 0,0-0,006) als in den bejagten (P = 8,9-15,6; He = 0,010-0,044) Popu¬ 
lationen. Die mtDNA (Cyt Z7)-Sequenzen zeigten denselben Trend. Diese Ergebnisse werden wie folgt 
interpretiert: Wahrscheinlich entwickelte sich der Fuchs im Laufe seiner Stammesgeschichte unter 
dem Druck von groBeren Predatoren, die im territorialen Fortpflanzungsgefiige einer Fuchspopula- 
tion eine gewisse Durchmischung bewirkten. In Abwesenheit dieser Rauber wurde derselbe Effekt 
moglicherweise durch die Jagd hervorgerufen und so eine groBere genetische Variation und damit 
Uberlebensfahigkeit bewahrt. Dies mag eine Ursache daftir sein, daB intensiver Jagddruck bisher 
nicht zu einer bedeutenden Reduktion von Fuchsbestanden gefiihrt hat. 
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